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Smooth estimation of mean residual life 
under random censoring 

Yogendra P. Chaubey* 1 and Arusharka Sen* 1 

Concordia University 

Abstract: We propose here a smooth estimator of the mean residual life 
function based on randomly censored data. This is derived by smoothing the 
product-limit estimator using the Chaubey-Scn technique (Chaubey and Sen 
(1998)). The resulting estimator does not suffer from boundary bias as is the 
case with standard kernel smoothing. The asymptotic properties of the estima- 
tor are investigated. We establish strong uniform consistency and asymptotic 
normality. This complements the work of Chaubey and Sen (1999) which con- 
sidered a similar estimation procedure in the case of complete data. It is seen 
that the properties are similar, though technically more difficult to prove, to 
those in the complete data case with appropriate modifications due to censor- 
ing. 

1. Introduction 

Let T denote a non-negative random variable representing the lifetime of a subject 
or a component in the context of survival or reliability studies. Given that the 
component has survived up to time t, the expected remaining life is given by the 
so-called mean residual life (MRL) function 

(1.1) m(t) = E(T t\T >t)= Jt s^ dU HS(t) > 0), 

where S(t) denotes the survival function (SF) of T. 

Yang [22] was the first to propose estimating m(t), by replacing the survival 
function in (1.1) with the empirical survival function (ESF). The reader may refer to 
Csorgo and Zitikis [6] and references therein for the vast literature on the properties 
of this estimator. There are several other articles that address the estimation of 
m(t) from a random sample of censored lifetimes. For example, Yang [21], Hall and 
Wellner [12], Ghorai et al. [9] and Gill [11], among others, studied the properties of 
estimators of m(t) based on randomly right-censored data, obtained generally by 
replacing S(t) by its nonparametric estimator such as the Kaplan-Meier estimator 
or suitable modifications thereof. However, due to the inherent discontinuity of 
the MRL estimator described above, there is a natural interest in deriving smooth 
estimators of the MRL function when the random variable T is assumed to admit 
a probability density function. In the complete-data case, Ruiz and Guillamon [17] 
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considered a kernel-smoothing-based estimator, and more recently, Abdous and 
Berred [1] have studied the properties of a local polynomial-based estimator of 
m(t). 

Chaubey and Sen [5] proposed an alternative MRL estimator in the complete- 
data case, based on a new smooth estimator of the survival function due to Chaubey 
and Sen [2]. Earlier, this smoothing technique was applied to hazard and cumulative 
hazard function estimation by Chaubey and Sen [3] in the complete-data case. Even 
though the method has been extended to density, survival, hazard and cumulative 
hazard function estimation for randomly censored data (Chaubey and Sen [4]), the 
MRL estimator has not been studied in the latter case. 

Chaubey and Sen [5] realized that the weighting scheme proposed for smooth 
estimation of survival and density function was not appropriate for the smooth 
estimation of the MRL function and modified these weights. The purpose of the 
present paper is to adapt the above method for smooth estimation of mean resid- 
ual life under the random censorship model. In this case we observe the censored 
lifetimes Zi = min(Tj,Cj) and the indicator variable Si = I(Ti < Cj), where C\, 
C*2 , . . • is an independent sequence of non-negative independent random variables 
representing random censoring times. The Kaplan-Meier product-limit estimator 
(PLE) for the survival function is defined by 

(1.2) Sn(*)=II !- ^ 



1=1 



i + 1 



where Z\-. n < ■ ■ ■ < Z n:n denote the ordered statistics corresponding to the obser- 
vations Zi, i < n and 5r, :n i is the value of S corresponding to Zi :n . 

For the censored data, kernel-based smoothing has been adapted for density 
estimation as well as for hazard function estimation. For example, the kernel-based 
smooth estimation of the density function is given by (cf., for example, Mielniczuk 
[15]) 

(1.3) /„(*) = K)- 1 r k (—) dF n (s), 



where F n (-) = 1 — S n (-), uj n (> 0), known as the band-width, is chosen so that 
lu h — > but nuj n -t oo, as n — > oo; fc(.) is termed the kernel function and is 
typically assumed to be a symmetric pdf with zero mean and unit variance. On 
the other hand, hazard smoothing is based on the Nelson- Aalen estimator of the 
cumulative hazard function H(x), namely, 

H n {x)= £ 



n — i + 1 

i:Z i:n <t 



as given by 



h n (t) = K)" 1 / k (—} dH n (s) 
Jo V w « / 

t-f n - i + 1 V w « / 

A list of available results concerning the above smooth estimator and other variants 
is contained in the survey article by Padgett and McNichols [16]. 
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Each of these methods can be used to derive a smooth estimator of mean resid- 
ual life. The smoothing technique considered here differs from the above standard 
techniques in such a way that it does not suffer from problems faced by kernel 
smoothing for non-negative data, and derives the estimator of the MRL function 
as a plug- in estimator. The technique presented here could also be used in general 
to smooth the Nelson-Aalen estimator, but this is not considered here. 

The form of the smooth estimator of the MRL function considered here is given 
in Section 2 and its asymptotic properties are investigated in Section 3. We also note 
that the results obtained here can be used to obtain the corresponding properties 
for the complete data case. In this process we have corrected the variance formula 
for the limiting asymptotic distribution of the resulting estimator. 



2. A smooth estimator of the MRL function 



As in Chaubcy and Sen [4], the following representation of the product limit esti- 
mator may be used for computation. Let M n be the total number of failure points 
and let the ordered failure points be denoted by Z* n , i = 1,2,..., M n , then we 
have 

(2.1) S n (t) = [] ( n ! fc fc ; i ) , Z* :n < t < Z* +Un , < j < M n 



th 



where h is defined by Z* n = Z kzin (Zq.„ = 0, Z\ u+l . n = oo), i.e. Z ki . n is the k, 
order statistic which corresponds to the i th failure. 

For smooth estimators of the survival function and associated functionals in the 
random censorship case, Chaubey and Sen [4] consider the following modification 
of S n (t) as in Efron [7]: 

(2 2) Sn ^ = iS n (t), for Z*. n <t< Z* +1:n , < j < M n 

I for t > Z n . n 

where for M n < n, we put Z^, In+1 . n = Z n - n . With this modification the smooth 
estimator of the survival function proposed by Chaubey and Sen [4] is given by 



(2.3) S n (t) -^5„(fc/A„) Pfc (A r; 

where 



k=0 



Pkifi) =e-"^j-,* = J l,2,... J 

X n being a constant to be chosen suitably. The seemingly infinite series in (2.3) is 
in fact a finite sum with the last term being the iV th term where N = [A„Z Tl:Tl ], [x] 
denoting the largest integer less than equal to x. 

Chaubey and Sen [4] argue that the data-adaptive choice of A n as given by 
A„ = n/Z n:n may be appropriate for estimating the survival function. For studying 
the asymptotic properties of the resulting smooth estimator of the MRL function we 
will make a deterministic choice of A„, which may be obtained by cross-validation 
(see above). 

Consistency of S n (t) can be proved using the same techniques as in Chaubcy and 
Sen [4]. Towards this we extract the following theorem and corollary from Stute 
and Wang [20]. 
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Let F and G denote the distribution functions of the survival time T and cen- 
soring time C and H denote the common distribution of Zi,i = 1,2, ...,n, i.e. 
1-H = Set 

t h = mf{x : H(x) = 1}, 

and write F{a} = .F(a) — F(a~). Finally let A denote the set of all atoms of H. 
Theorem 2.1 (Stute and Wang [20]). Let F n = 1 — S n . Suppose that F and G 
have no jumps in common, and <p is F-integrable, then we have with probability one 
and in the mean, 

(2.4) lim ( 4>{x)F n {dx) = [ cj>{x)F{dx) + l {TH&A) <l){T H )F{T H } 

J J x<th 

(where the second term should be zero if Th ~ oo). 

Two special cases of this result are of importance for further analysis and are 
given as corollaries below. 

Corollary 2.1. Suppose F and G have no jumps in common, then we have 
sup \S n (t) — S*{t)\ — > with probability 1, 

t<T H 



S*(t) 



where 

S(t) if t < T H , 

S(t h ) + 1 Th <= a F{t h }, if t>r H . 

The uniform consistency of S n holds on (0, th) if either F{th} = or F{th } > 
and G{t h } < 1. 

Remark 2.1. We will consider both F and G to be absolutely continuous hence, 
the uniform consistency for the product limit estimator will be assumed to hold on 
the support of H. In addition, if th = Tp and F is continuous (see Remark 3 of 
Stute and Wang [20]) then we have 



n — >oc 



(2.5) lim / <j)(x)F n (dx) = / cj>(x)F(dx) 



This implies the following corollary. 
Corollary 2.2. If th = Tp and F is continuous with J \x\F(dx) < oo, then with 
probability one we have 

(2.6) lim / (x-t)F n (dx)= (x-t)F(dx), 

n ^°°Jx>t Jx>t 

where the above convergence is uniform on the interval [0,r] for any r < th- 
Remark 2.2. If an addition F is concentrated on non-negative values, the above 
implies that with probability one, \im n ^oo J x>t S n (dx) = J x>t S(dx) uniformly on 
the interval [0, r] for any r < tjj- 

Now we are ready to establish the strong convergence of S n (t). 
Proposition 2.1. Let H and F have the same support and assume that F is 
absolutely continuous. Consider a sequence of positive numbers A„ — ► oo as n — > oo 
then we have 

sup \S n (t) — S(t)\ — > 0, as n — ► oo. 

t<T H 

Proof. The estimator considered here is based on the following theorem (see Lem- 
ma 1 in Feller [8], pp. 219): 



Smooth estimation of mean residual life 



39 



Let u(t) be a bounded and continuous function on R + . Then 

(2.7) e~ At u{k/X)(\t) k /k\ -> u(t), as A ^ oo, 

fe>0 

uniformly in any finite interval J contained in R + . 

Using the above result, since S(t) is bounded, continuous, nonnegative and non- 
increasing (on R + ), we have 

(2.8) S*(t) = e~ Xnt S(k/X n )(Xt) k /kl -► as A„ -► oo, 

fc>0 

uniformly in any finite interval J C R + . 
Now we have 



sup {\S n (t)-S* n (t)\} < max |5„(fc/A„) - S(k/\ n \) 

t<T H {k:{k/X n )<T H } 

(2.9) +5„(r H )Prob{iV„ > [r H A„]}, 

where iV n is a Poisson random variable with mean t\ n . Both terms converge to zero 
with probability one, the first one due to Corollary 2.1 and the second one because 
Sti{th) is bounded and Prob{iV rl > [TffA ra ]} — > as A n — > oo. For analyzing the 
difference between S n (t) and S(i), we see that 

sup \S n (t) - S(t)\ < sup \S„(t) - S*(t)\ + sup \S*(t) - S(t)\. 

t<TH t<TH t<TH 

Using (2.8) and (2.9) we find that both terms on the right-hand side of the above 
inequality converge to zero and hence the theorem holds for th < oo. For the case 
of th = oo, the second term on the right-hand side of (2.9) is exactly zero and the 
theorem follows similarly. □ 

The estimator of the MRL function based on the Kaplan-Meier estimator as 
given by Yang [21] can be written as 

r°° g (x)dx 
m n (t) - Jt "V I(t < Z Mn:n ). 

We can similarly define the smooth estimator of the MRL function as 

f x >t( x ~ t)Fn(dx) J t °° S n {u)du 



(2.10) rh n {t) 



S n (t) S n (t) 



where F n = 1 — S n . 

The computational form of the above estimator is now provided in the following 
proposition as given in Chaubey and Sen [4]. 

Proposition 2.2. The smooth estimator of the mean residual life as given by (2.10) 
can be represented as 



(2.11) m„(t) = 



1 l^k=0 l^r=0 (fc-r)! X„ ' 

K T N <* A -) fc s ( k ) ; 
2^k=o k\ °»IaJ 



where N = [X n Z n:n ]. 

The representation given in (2.11) follows by direct integration and the fact that 
S n (x) = for x > Z n:n . 
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3. Asymptotic properties of rh n (-) 

In this section we present the asymptotic properties of rh n (-). We establish its strong 
uniform consistency in Theorem 3.1 and asymptotic normality in Theorem 3.3. 
Theorem 3.1. Let F and G be absolutely continuous with common support and 
no common jump points such that m(t) < oo for all t € R + . Then for any compact 
interval C C 1Z + , we have for X n — > oo as n — ► oo 

(3.1) || rh n — m \\c= sup \rh n (t) — m(t)\ — > a.s. as n — > oo. 

tec 

Proof. If F is compactly supported, i.e. S(t) = for all t > to(< oo), then m(t) 
is not defined beyond to and for all t < to, it is bounded. Thus we find that m(t) 
is a Hadamard diffcrentiablc functional of S(t). Hence, we can claim the same 
convergence property for m n (t), as that of S n (t) due to Proposition 2.1. Therefore, 
we confine ourselves to the case of infinite support, i.e. t# = oo. 
First note that (2.11) can be written as 



(3.2) rhn(t) 



S n (t) 



where Pk(x) = ^j<j.Pk{x)- Writing the numerator in the above equation as 
Y,k G n( k /X n )p k (tX n ), where G n (k/X n ) = (1/A„) J2j<k S n{j/K), we note that 

cJ±\-± s J±\<c°(±-)<a/ k 



A n / A ri \ A n J \ X n J \ X„ 

where 

S n (x)dx. 

Hence, 

(3.4) sup \G° n [A J - G n (^-] | -► a.s. asn->oo, 

where C = [0,r],r < oo. Further, since, G„(t) is non-decreasing and bounded a.s., 
using Hille's theorem, we get 

(3.5) ^2pk{tX n )G n {k/X n ) -> G°(t) a.s. as n -> oo Vi £ C. 

fc>0 

Now using Corollary 2.2 (see (2.6)) we have 

/•OO 

(3.6) G„(t) — > / S{x)dx — m(t) a.s. asm oo. 



Combining (3.5), (3.6) along with the fact that S n (t) 1 — ► 5(i) 1 , uniformly in 
any compact interval we get the result in (3.1). □ 

Before establishing the asymptotic normality of rh n (t) we need the following 
notation. Let 

H Q (z) = P(Z <z,S = 0)= f Z S(y)G(dy), 

Jo 
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and 

H x {z) =V(Z<z,8=l)= f (1 - G(y))F(dy). 

Jo 

Define 

7o(x) = cxp Uo T^)l' 
i r 00 

71 (x) = 70 (x) = - _ ^ y I(x < w)(/){w)'yo(w)H 1 (dw) 

and 

/ ^ _ r / s f°° f°° I(v < x,v < w)(/>(w)y (w) . . , . 
72 (a;) = r (x) = y ^ (l-g(t;)) 2 H (dv)H 1 (dw). 

Under the following assumptions, 

(3.7) J 4: 2 (x)^(x)H 1 (dx) <cx) 
and 

(3.8) I |0(.T)|C* 1/2 (x)F*(dx) < oo, 
where 



C(x) 



G(di/) 



/ (1-G(j/))(1-H(y))' 

Stute [19] proved the following CLT. 
Theorem 3.2. Under (3.7) and (3.8) 

J cf>d(S n -S*)^N(0,al), 

where 

g% = Var{0(Z) 7o (Z)<5 + 7l (Z)(l -S)- l2 (Z)}. 

We show here that a similar central limit theorem holds for the smooth estimator 
rh n (t). As mentioned in Remark 2.1, we suppose that G is also continuous and 
assume the conditions (3.7) and (3.8) with 4>(x) = x, i.e. 

(A.l) 



7 1 - G(x) 



and 

(A.2) 



DC'. 



where C{x) now becomes 

G(dy) 



G(x) 



/ (1-F(y))(l-G(j/)) 
Note that the above two conditions imply 
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(B.l) A.l holds, replacing x 2 in A.l by (4>£{x)) 2 or (</>f (x)) 2 and A. 2 holds re- 
placing (the first) x in A. 2 by <ftf(x) or 4>f{x) for any < t < th, where 
<tf{x) = (x- t)I(t > x), $(x) = I(x > t). 

(B.2) Given e > 0, there exists M, t < M < th, such that 

TH X !^^ M) F{dx) < e and / " xl(x > M)y/C{x)F(dx) < e. 
1 - G{x) Jo 

The above equations are similar to (2.6) and (2.7) of Stute [19] (see p. 435). 
Further, note that 
(B.3) for any < a < t h , A.l and A. 2 hold for 4>{x)I{x < a), i.e. 



f™ <j) 2 (x)I(x <a) 
Jo 1 - G(x) 



and 



\<j){x)\I{x < a) v / C(x)F{dx) < oo, 

Jo 

whenever / cj) 2 {x)F{dx) < oo. (sec (2.3), p. 432 of Stute [19]). 
Theorem 3.3. Under assumptions A.l and A. 2 and if X n — > oo, for any t < th 

Vn~{m n (t) - m{t)) ► JV(0, cr 2 (t)), 

where 

a it) = asymptotic variance of {\/n(m n (t) — m(t))}. 

(See Yang [21] or (3.11) below). 

For further analysis we will use the following representation of the K-M estimator 
of the survival function as given in Stute [19]: 

n 

s n (t) = Y,w m i(z l:n >t), 

i=l 

where for 1 < i < n, 



W in = 



6, 



i:n] 



i— 1 

(n - * + i) n 

v .7=1 



n-j 



i s 



[j:n] 



is the mass attached to the ith order statistic Zi- n under F n . Therefore the smooth 
estimator of the survival function may be written as 



where 



S n (t) =J2W in <t>% t (Z i:n ), 



°° ( k \ 

(x) = ^ W (U„)I 

fc=0 



P(T 1 <»). 

An, 
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13 



S„(y)dy 



where 



i=l 
n 

i=l 



V J(Z i:n > — ) / Pk(Ky)dy 
k=o An Jt 

00 fc f 1 fc 

5^ I(Z iin > — W — ^>,(^ 



,fe=0 
00 



r=0 



E( [A " Zt: "l fc+1 ) 



fc=0 



1 A 1 \ ^[A n x]-fc + l / fc 

= 2^ 1 4 Pk{t\n)I [ X > 



k=0 



An 

[A n x] - iV„ + 1 



A,, 



I(N n < xX n ) 



First note that, as n — > 00, 



<^ t (x) -> J(t < a) + -I(t = x) and 4>i t {x) -> (x - t)I(t < a) 



Let us therefore consider, 



(3.9) 
(3.10) 



S„(t) = X^ m 0f(Z l: „) 

i=l 

/ S n (y)dy = y;W in ^(Z i: „), 



where 



<£f (x) = /(x > i) and (j)f(x) = (x - i)/(x > i). 
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We have by Theorem 3.2 

\fn{m n (t) - m(t)) 



' \ S n (t) 

- (IT S n (y)dy - / f °° S(y)dy) 
S n (t) 



^ Sjt) 



— y/nm(t) 

A 



(S n (t)-S(t)) 

S n (t) 



— \/nm(t) 
1 



S n (t) 



S(t) 



-E F (cf>? 



- G(Z. t ) 

n 

-^E(( 1 -^( z *)- r ^( z *) 

* i—l 

\fn V 1 



G(Zi 



m(t) 



(3.11) 



V i—l 

o P (l). 



The claim in the above equation is justified by the assumptions in (3.7) and (3. 
In view of the above, to prove Theorem 3.2 it is enough to show that 



(3.12) R n (t) = 



J7° Sn(y)dy / t °° S n {y)dy 



S n (t) S n (t) 
To establish the above equation, we write 

Vn[Jt°° S„ {y)dy - J t °° S n (y)dy 



op(l) as n — > oo. 



Rn{t) 



S n (t) 

■ [SrM - SrM ( fT Sn(y)dy \ 



S n (t) 



S n (t) I 



(3.13) 



R 



*M ( St 00 S n (y)dy \ R*(t) 

S n (t) 1 S n (t) J § n (t) 



say. 



To complete the proof of the assertion in (3.12), we establish the following lem- 
mas. 

Lemma 3.1. Under A.l and A. 2, 

Rn(t) = op (I) as n — > oo. 
Lemma 3.2. Under A.l and A. 2, 

Rn W = °p(1) as n ~* 00 ■ 
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Proof of Lemma 3.1. We have 



Rn(t) 



\M/ S n (y)dy- \ S n {y)dy] 



°° f°° ( k \ 

^2 Vn J [S n y—j ~ S n(y)]Pk(Ky)dy 



fe=0 



E 



fc=0 



a n y— | - o> n (y)]pk(\ n y)dy 



00 />oc / k \ 

V^E / t 5 hr -S(2/)H(A„y)^, 



where 

We may therefore write 



M = MSn(y) - S(y)). 



E 

fc=0 



E 



M 



at n {y)]pk{Ky)dy 



a„ ( — ) - a„(y) 



Pk(Ky)dy 



E 

fc=0 



S(-)-S(y) 



Vk{Ky)dy 



(3.14) 



where M > M\ and i < A/i < t# is chosen to satisfy B.2 for a given e > and 
A n (t) is a deterministic function converging to zero, as shown in Chaubey and Sen 
[4]. Next, we have 



r n ,M(t) = 



A I 



E 



k:\ -g--y\<h„ 

M 



&n(y) 



E 



k:\-fi--y\>h n 



a n [ — )- a n (y) 



Pk{\ n y)dy 
Pk(Ky)dy 



A I 



i(y)dy 



A I 



r n.M.2(y)dy, say. 



By Corollary 3.1, p. 1316 of Deheuvels and Einmahl [10], 



(3.15) 



SU P \ r nMi(y)\ = 0(h n [- log h n + log log n]) 

t<y<M 



provided h n I 0, ^ -»• 00 and -> 00. Further, by (14), Inequality 1, p. 485 

of Shorack and Wellner [18], we have 



(3.16) 



N 

\ M ,2(y) ^ v^p[|-r-tfl >k 

An 

A rl h n 



< v n exp I — 
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a.s. 



It follows from (3.15) and (3.16) that 

< M (i) < (M -t)0(h n [- logh n + loglogn}) 

+ (M - t)y/n cxp 

(3.17) — ► a.s. as n — > oo 

for any suitable choice of /i„. 

Next consider i?^ M (i). We have 



2M 



00 /"oo / k \ 

E / I a « ( I - ) ~ ot n {y)]p k {X n y)dy 

k=0 \ An ' 



k=0 



\/™E j M [ S n y—j - S y—j}pk(X n y)dy 



OO »x 

<E / [<Sn(y) ~ S{y)]p k {X n y)dy. 

k=0 JM 



Now using the fact that 



/ Pk(X n y)dy = -Vp r (A„M) 



in the first term on the RHS of the above equation, we have 

([X n Z i:n ] + 1 - r) 



^n,Af(*) 



E 



E^» 



?:=i 

([A„X] + 1 - r) 

An 



An 



> T-) 

An 



/(X > — ) 

An 



p r (A n M) 



E WinCim ~ M)I(Z i:n > M) - E F {{X - M)I{X > M)) 



i=l 

* E t 

r:(r/A„)<Mi 



* E i 

r:(r/A„)>Mi 



(3.18) 



-y/nj <t>t I {x){F n {dx)-F{dx)) 



where i < Mi < A/ < th are as given after (3.14). 

Note that the summands in M 1 and i?„ m 2 are of the form 



where 



n / ip r>n {dF n - dF) Pr (X n M), r = 0, 1, 2, . . 



Wr,n\x) = " I(X > (r/Xn)) 

An 



- — I (x > (r/A„)) + (x > (r/A„)) , 



JV 



X r , 



Smooth estimation of mean residual life 



17 



where < a„(x) < 1. Furthermore, we have 

i>r,n{x) = (Pn.(r/\ n ){x) + — 1 {x > (r / \ n )) 

A n 

(3.19) < xl(x > Mi) + {l/\ n )I(x > Mi) for (r/A„) > Mj 

(3.20) < x + 1 for all r = 0,1,2,.. . 

Using the bounds in (3.19) and (3.20), and arguing as in the proof of Theorem 1.1, 
pp. 435-437 of Stute [19], we get by B.2 that 



(3.21) y/n J ' ip rtn (dF n - dF) = O p {\) uniformly for < r < A n Mi 

(3.22) ^J^, M F„- iF) . Op{e) uniformly for r > A„Mi 
(dF n -dF) = P (e) 



(3.23) 
Thus 



>n I <p 



MA 



VTl[J ^r,n(dFn-dF)}p r (\ n M) 



r:(r/\ n )<M- 



P (l)Prob{^^- < Mi}, where N n , M ~ Poisson(A„M) 

An 



(3.24) 



P (l)Prob 



N n , M - A„M r—Ah - M 

S V 



M 



as n — > oo, 



since Mi — M < and A„ — » oo as n — > oo. 
Also 



-ft;. 



M.2 



Vn[ J AAdFn - dF)]p r (\ n M) 



(3.25) 



r:(r/i„)>Mi 

Opto. 



and Ri >M> z = P (e) by (3.23). 

Using (3.17), (3.18), (3.21)-(3.23), (3.24) and (3.25), Lemma 3.1 is established. 



□ 



Proof of Lemma 3. 2. 



E 

oo 

E 

k=0 



S n (—)-S n (t) 

An 



a n(i—) ~ at n (t) 

A n 



E 

fc=0 



S(-)-S(t) 

An 



Pk{A n t) 
Pk(A n t) 

Pk(A n t), 



so that negligibility of R^(t) follows by arguments similar to, but much easier than, 
those in the proof of Lemma 3.1. □ 
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Now by (3.13), R n (t) is asymptotically negligible and hence the asymptotic nor- 
mality of ^/n[rh n (t) — m(t)] is now established and this also shows that the limiting 
distribution is the same as that of y/n[m n (t) — m(t).] 

Remark 3.1. In view of the above asymptotic result we note that a similar result 
holds in the complete data case, where 5=1, G(-) = 0, r^(-) = 0, so that by 
(3.11) the expression for the asymptotic variance simplifies to 

a 2 {t) 

= var ([(X - t)I(X > i) - m(t)I{X > t)]/S(t)) 
= (5 2 (t)) _1 [var((X-t)/(X > t)) + m 2 (t)v&r(I(X > t)) 
-2m(£)cov((Y - t)I(X > t), I(X > t))} 

= {S 2 (t))- l [E{(X-tfl{X>t))-{ S{y)dy) 2 +m 2 {t)S{t){l-S{t)) 



-2m(f)(l - S(t)) I S(y)dy] 

{S 2 ^))- 1 [E((X - t) 2 I{X > t)) - m 2 (t)S 2 (t) + m 2 (t)S(t)(l - S(t)) 
-2m 2 {t)S{t){l- S{t))l since / S(y)dy = m(t)S(t), 



= (S 2 (t)rM£((* - tfl{X > t)) - m 2 (t)S(t)} 
= var((X -t)\X >t)/S(t), 

which is the same as the one obtained by Yang [22]. Hence, for the complete data 
case we have 

V7i{rh n (t) - m(t)) -> N(0, a 2 (t)) as n — > oo. 
This corrects the asymptotic distribution obtained in Chaubey and Sen [5]. 

4. Concluding remarks 

We propose a smooth estimator of the mean residual life function under random 
censoring and investigate its asymptotic properties. To the best of our knowledge, 
no such estimator has been studied so far in the random-censoring case. We adapt 
the smoothing technique of Chaubey and Sen [2, 4], which is appropriate for non- 
negative data, especially in terms of avoiding boundary bias. We establish strong 
uniform consistency and an asymptotic weak representation of our estimator as an 
average of iid random variables. The latter, which automatically yields asymptotic 
normality, is established using the results of Stute [19] and Deheuvels and Einmahl 
[10]. The asymptotic variance of the estimator for uncensored data (Chaubey and 
Sen [5]) is derived as a special case. 

Acknowledgments. The authors would like to thank a referee and Mervyn J. 
Silvapulle for their comments which have improved the readability of the original 
manuscript. 
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